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Abstract 

The existence and uniqueness of a solution to a generalized Blasius equation 
with asymptotic boundary conditions are proved. A new numerical approximation 
method is proposed. 
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1 Introduction 

We study the BVP of the form: 

x'" + cx p ■ x" = 0, x (0) = = x'(0), \imx'(t)=(3, (1.1) 

t— >oo 

where p > 1, c and are positive constants. The problem is motivated by the classical 
Blasius equation describing the velocity profile of the fluid in the boundary layer where 
c=|,p = /3 = l. The Blasius equation is a basic equation in fluid mechanics which 
appears in the study of the flow of an incompressible viscous fluid over semi infinite 
plane. Blasius ([!]) used a similarity transform technique to convert the partial differential 
equation into his famous ordinary differential equation 

x'" + \x ■ x" = 0, x{0) = = x'(0), ]ima/(t) = l, (1.2) 

t— >oo 

where x is the stream function x = , U is the fluid velocity, v is the fluid kinematic 

viscosity and t is the similarity variable defined as t — y<i\J where y%, yi are Cartesian 
coordinates with y\ pointing along the free stream direction and y<i perpendicular to y\. 
We refer to [2J [3] for an excellent introduction to the problem. A series expansions 
method was used to solve (11.2|) by Blasius. There has been appeared many analytical and 
numerical methods handling this problem since the Blasius's work, [61 [8] for instance. 
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In the first part of this paper, the existence and uniqueness of (11. ip will be analytically 
proved by changing the boundary value problem to an initial problem. Using the obtained 
estimates we will be able to find the value of a = x"(0) which guarantees that the solution 
x a on an initial problem 

x" + cx p ■ x" = 0, x(0) = = x'(0), x"(0) = a, (1.3) 

is the solution of (11. ip we are looking for. In the second part of the article, a new numerical 
approximation method is proposed. 



2 Auxiliary lemmas 

Let x a stand for the unique solution satisfying initial conditions 

x(0) = = x'(0), x"(0) = a. 

If a < 0, then x a and x" a are negative for small £'s thus the solution is concave and negative 
for all arguments and it cannot solve (jl.ip . For a = we have a trivial solution x a = 
and the seeking solution can be obtained for a > 0. 

Lemma 1. TTie x a is defined for all t > 0. 

Proof. x" a cannot vanish at any point t by the uniqueness of solutions of initial 
value problems: x(t) = C\t + solves our ODE. Hence dividing the equation by x" a and 
integrating on [0, t] we have 

x" a (t) = aexp ^— c J x a (s) p ds^j , (2.1) 

which implies 

x ' a (t) = a exp x a (T) p dT^j ds. (2.2) 

Integrating once more and applying the Fubini Theorem we get 

-t 



x a (t) = a J (t — s) exp ^— c ^ x a (r) p c?r^ ds. 



(2.3) 



By (I2.ip . (12.2p and ( 12.31) we have apriori estimates: 

< x a (t) < -at 2 , < x'Jt) < at, < x"(t) < a 
for any t > 0. It follows ([7], p. 146) that x a is extendable to [0, oo). 

Lemma 2. For any a > 0, lim^oo (£) = and there exists a finite and positive limit 
h{a) := finite x' a (t). 
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Proof. Since x" a > and x' a > 0, then \im t ^ +oc x a (t) = +00. Moreover, since x" a > 0, 
x a > and = —cx p a x'^ then x" a is a decreasing function so lim^oo x"(t) = g a G [0, a). 
Suppose g a > 0. From x' a " = —cx^x'^ we get that 

lim <'(*) = -00. (2.4) 

t— >+oo 

On the other hand 

V* > 3 s t 6 (i, t + 1) x' a '(t + 1) - = <'(s t ) 

and hence lim 4 _;. +OC) x"'(t) = 0, which contradicts (12.41) . 

Since x" a > 0, then x' a is an increasing function (so the limit defining h(a) exists, 
possibly infinite). From \im t ^ +oc x a (t) = +00 we get there exists t a > such that 
cx a [t) p > 1 for t > t a . From (12. 2p . we obtain for t > t a , 

x' a (t) < a J exp ^— c J x a (r) p dr^j exp — J dr^j ds < 

rta ft POO 

< a ds + a e " (s "* a) ds < at a + ae ta / e~ s ds = a(t a + 1). 

JO Jt a Jt a 

Thus 

h(a) < a(t a + 1) (2.5) 

for any a > 0. 

Lemma 3. For any a > 0, there exists a finite and positive limit 

fi(a) := limt_ >0O (/i(a)t — x a (t)). It means that the graph of x a has a slant asymptote and 
the following estimates hold: 

max(0, h(a)t — //(a)) < x a (t) < h(a)t. (2.6) 

Proof. The function t h-> h(a)t — x a (t) is increasing, hence the limit from the assertion 
exists but it can be infinite. Suppose it equals +00. By the arguments from the proof of 
the lemma[2j we have x'"(t) < —x'^(t) for t > t a . Integrating this inequality from s to +00 
and using the fact x' a '(+oo) = 0, we get 

-x" a {s)<-h{a)+x' a {s). 

Next integration from t a to t leads to the following inequality 

x' a (t a ) - x' a (t) < -h(a)(t - t a ) + x a (t) - x a (t a ) 

or equivalent ly 

h(a)t - x a {t) < h(a)t a - x a (t a ) + h(a) - x' a (t a ). 

Thus 

< ri(a) < h(a)t a - x a (t a ) + h(a) - x' a (t a ). (2.7) 
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The last part of the assertion is a simple consequence. 

For an upper bound on h(a), //(a) depending explicitily on a we use x a (t) < at 2 /2 to 
flZS]). Hence, 

x a {t) >aj^(t-s) exp (- J c^r 2p dr^j ds = a j^t-s) exp (-c ^ 1)2? / P+1 ) da. 
One can easily show that the function 

¥>(*) : = y (* - s ) ex P ds, fc = 2P ( 2p+1 y « = 2p + 1 

has a similar behaviour as x a in the sense that its graph has an asymptote x = h*t — jj* , 
where 

h* = exp (-ks a ) ds = ^M, = jf s exp ds = (2.8) 

and its graph sits above this line. Hence, 

x a (r) > a 1 -^ • r(^ I )(^ I?? ) W . r _ a i-^i . r(^ T )(2p + 1) W(f 

p+i i 
^o(t) > c 2 ■ a 2 f +1 • t - c 3 • a 2 p+! , (2.9) 

where 

^==r(^i)(^i^)^, c 3 :=r(^ I )(2p + l)^(f)^i. (2.10) 
Now, we are able to get appropriate estimates for h(a). 
Lemma 4. For any a > 0, 

c 2 ■ a ( p+1 )/( 2p+1 ) < h(a) < Cl ■ a ^/^ +1 \ (2.11) 

where 

. _ c 3 | r(i/(p + i)) 

1 ' C 2 cVCP+l) . ( C2 (p + I))p/(P+1)- 

Proof. For a lower bound we apply the estimate x a (r) < |ot 2 to the equality 

/i(a) = a J exp ^— c ^ x a (r) p dr^ ds. (2-12) 
This leads to the inequality 
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(12. 8p for = 2P (2p+i) > a = 2p + 1 gives the lower bound on /i. 

For an upper bound on h we use the lower estimate of x a - (12.91) to the equality (12.121) 
and we get 

pc 3 /c2-a~P/ < - 2 P+ 1 1 

h(a) < a ds+ 



f°° ( r f p+i i v \ 

/ exp I — c / ( C2d 2p+1 T — c s a 2 p+ 1 ) dr ds 

ic 3 /c 2 -a-P/( 2 P+ 1 ) V Jc 3 /c 2 -a-P/( 2 P+ 1 ) ^ ' J 

C3 p+i 1 p Z" 00 / C P+i „,i\ , 

< — a 2 ^ 1 H a 2 ^ 1 / exp ( -- — — - — a 2 p+H p+ ) dt, 



p+i 



where we used linear substitutions twice. At last, using (I2.8P for k = j^p^a 2 p+ 1 
a = p + 1 we have for any a > 0, 

h(a) < ( c ± + r(i/(p + i)) \ ^ 

W S U CV(P+D ■ ( C2 (p + 1)) P /(P+D ) 

The next lemma presents estimates for n{a). 

Lemma 5. For any a > 0, constant fi(a) satisfies the following estimates: 

c^-a 1 /^ <^a)<c 5 -a^ 2p+1 \ 

where 

2 2p/(2p+i)r( 2 /(2p+l)) 



c 4 := 



(2p + l)(2p-l)/(2p+l) c 2/(2p+l) 



Proof. From (Q and (l2TT2|) we get 

/i(a) = a sexp I — c / x a (r) p dr I <is. (2-13) 



Using the estimate x a (r) < ^ar we get 



2 

cu 2p+l 



a(a) > a s ■ exp ; -s p ds. 

^ K ' - Jo v 2f(2 P +i) ; 

Using ( 12. 8 p for k = 2 p(2p+i) ' a = 2p + 1 we obtain the lower bound. On the other hand, 
from (El), d2H3D we have 

/•c 3 /c2-a-P/( 2 P+ 1 ' 

/i(a) < a sds+ 



f°° ( r / p+l \ 

1 / sexp I — c / ( C2a 2 p +1 r — c 3 a 2 p +1 ) dr J ds < 

Jc 3 /c2-a-P/( 2 P+ 1 ) ' V Jc 3 /c 2 -a-P/( 2 P+ 1 ) ^ ' ) 

C\ -J- 1 -zL. /"°° / C _ P±i. , 

< -4ra 2 p+i + ^a 2 p+i / texp -- r — a ^p+H p+ dt+ 

~ 4 J * V (P + 1)C2 / 



+ ^ / exp - - C . a-Wt^ +1 ) dt, 
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_ p+1 

where we used linear substitutions twice. At last, using (I2.8P for k = ( p+ i) C2 a 2p+1 , 
a = p + 1 we have for any a > we get the upper bound. 

3 Main results 

Now, we are able to prove the existence of a solution to (II. ip . 

Theorem 1. The BVP U.l\) has a solution for any (3 > 0. 

Proof. The function h : [0, oo) — > R is continuous on (0, +oo) by the continuous 
dependence of solutions of ODEs on initial conditions and locally uniform convergence of 
the integral 

J exp ^— c J x a (r) p dr^j ds. 
By the estimates from lemma HI we have 

lim h(a) = 0, lim h(a) = +oo. 

Thus, for any (3 > 0, there exists a > such that h(a) = (3. For (3 = 0, it is obvious. 

Finally, the uniqueness of the solution of (11. ip will be proved by using the ideas 
from [1]. For any a > consider the one-to-one function v a ■ [0, h(a) 2 ) —> [0, oo) such 
that v a {x' a (t) 2 ) = x a (t) for each t > 0. It is well defined since x a and x' a are increasing 
functions and it belongs to C 2 (0, h(a) 2 ). Substituting y = x'(t) 2 , we shall find an ODE 
satisfied by v (we omit subscript a for simplicity). 

x (t) = v(y), x'(t) = v'{y)2x'(t)x"(t) 

hence, 



x"(t) = — x"'(t) = - v " {y) ■ 2x'(t)x"(t) = - v " iy) ^y 

X[t) 2v'(yY X[ ) 2v'(y) 2 AX[t)X [t) 2v'{yf 
Put x and its derivatives in our ODE and find 

... , cv(y) p v'(y) 2 
v"(y) = yy, yy> . 3.1 

From boundary conditions on x we get 

v(0) = 0, v'(0) = —, lim v(y) = +oo. (3.2) 

2a y^h(a) 2 - 



Now, we are in position to prove 

Theorem 2. The solution of U.l\) is unique. 
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Proof. We need to show that the function h is one-to-one. Suppose that h{a\) = h(a 2 ), 
a 2 > a\ and take V\ and v 2 obtained by x ai and x a2 , respectively, that is Vi satisfies (13.11) 
with boundary conditions ( 13. 2p (for a = i = 1,2). Put w = v x — v 2 . Then w(0) = 0, 
w'(0) = > 0. Notice that w' > on the whole interval (0, /i(ai) 2 ) - both function 

are defined on the same interval. 

In fact, if it is not true, then there exists s in this interval such that w' > on (0, s) 
and w'(s) = 0. Hence w(s) > w(0) = and 

„. . , w'(s — £ ) — w'(s) 
w"(s) = lim — i ^- ^ < 0. 

On the other hand, 

w"(s) = v'l(s) - v' 2 \s) = cs-^Msf - v 2 {sY)v' ai {sf 

since w'(s) = implies v' a (s) = v'(s). But Vi(s) p > v 2 (s) p from w(s) > and this gives 
w"(s) > - a contradiction. Thus, we have w > and w' > on (0, h(ai) 2 ). 
Set Vi = l/v[, i = 1, 2 and W = Vi - V 2 . We have, for any y G (0, %i) 2 ), 

r M = - vM = _J*M + = c ^-oW < o 

My) w) Vv 



from w(y) > 0. Hence, 



and 



lim < 2(ai - a 2 ) < 0. 

^Ti(ai) 2 - 



On the other hand, Vi(x' a .{t) 2 ) = 2x".(t) implies 



lim Vt(y) = lim V^«(t) 2 ) = 2 lim <(t) = 



and, therefore, 
which contradicts the previous inequality. 



lim W(y) = 



4 Numerical approach 

All numerical methods cannot work on the infinite interval [0, oo) and we do not know 
the exact value of a = x"(0) for the solution. Our earlier results make possible to find a 
finite interval [0, T] for any positive value e of the error control tolerance such that 

x"(T) < e, h(a) - x\T) < e, x(T) - (h(a)T - //(a)) < e. 
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Since all these functions decrease, all three inequalities hold for any t > T. First, by using 
estimates ( 12. lip , we can find an interval [a m i n ,a max ] such that h(a m i n ) < j3 < h(a max ). 
Next, by (12.11) . we need 

a max exp (~ c J o x a min {. r f d ^ < e > 

by (12.121) . we should have 

a ma x J exp J x a mm (jY dr^j ds < e, 
and by(T23D and (12~T3|) . we get 

-a max j (T - s) exp f-c J x amm (r) p dr^j ds < e. 
We do not know the function x amin but we can use estimate ( 12.9ft to get 

<W exp (-c [ (c 2 ■ • r - c 3 • a^ 1 \ dr) < e, (4.1) 







exp -c / c 2 • a^ 1 • r - c 3 • a^ 1 dr ) ds < e, (4.2) 



/■oo { [ S f - E±1 - -^—\ P \ 

- a max (T - s) exp f-c jf ( c 2 • a^ 1 • r - c 3 • a^ 1 J dr J rfs < e. (4.3) 

We start with a family of initial value problems 

x a (t) for t G [0, T], x a (0) = = x' a (0), x'^(0) = a, a e [a min , a max \. (4.4) 

If we approximate this solution in [0, T] with an error less than e, then the best approxi- 
mation of x in [T, oo) is 

x(t) = 0t + (x a (T) - 0T) . 

The lower and upper bounds for the second derivative describe the shooting window - for 
each a the only solution in this direction exists at t = T, and the computed value of x'(T) 
is more and more close to the expected limit value /3. As long as /3 is contained between 
the computed values x'(T) of the best two shots, we apply the classical bisection method: 

If y and z are solutions such that y"(0) < z"(0), and there is y'{T) < (3 < z'(T), then 
the next problem to solve is (l4~4l) with a = (y"(0) + z"(0))/2. 

Examples. 

While solving the initial value problems we apply an adaptive Runge-Kutta-Fehlberg 
method RK45 [5], in which a tolerance parameter e controls local error of the method. 
The values of e range from 1CT 8 to 10~ 14 . For representing real values we use standard 
16 — 17-digits double data type. 
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Numerical results for the classical Blasius equation p = 1, c = 1/2, and [3 = 1. 

Here a min = 0.2694860459, a max = 0.3420953216. For T = 14 we get the all three 
inequalities (14. ip . (}4.2j) and ( 14. 3 p for e = l.Oe — 14. Below N stands for a number of steps 
in RK45 (average): 



e 


N 


a 


\x"(0) - a\ 


\l-x'(T)\ 


x(T) 


1.0e-08 


112 


0.332057330068201 


6.15e -09 


1.2e-08 


12.279212180321 


1.0e-09 


197 


0.332057335646357 


5.69e -10 


l.le-09 


12.279212327474 


1.0e-10 


338 


0.332057336149237 


6.60e -11 


1.3e-10 


12.279212340740 


1.0e-ll 


611 


0.332057336210248 


4.95e -12 


9.9e-12 


12.279212342350 


l.Oe- 12 


971 


0.332057336214903 


2.94e -13 


5.7e-13 


12.279212342472 


l.Oe- 13 


1831 


0.332057336215154 


4.19e-14 


6.5e-14 


12.279212342479 


1.0e-14 


3346 


0.332057336215186 


1.06e-14 


5.5e-16 


12.279212342480 



The last value of a differs in two last digits from the one cited in [2]: 



a = 0.33205733621519630. 

The last two columns of the table have been computed for e = 10~ 14 . As there is x"(T) = 
7.68e — 13, for t > T the straight line approximation of the solution is the most effective. 

Numerical results for the equation with p = 7, c = 1/2 and (3 = 1. 

Here, a min = 0.3733978388, a max = 0.3805482427. As above for the tolerance e = 
l.Oe — 14, the interval [0,4] is sufficiently large and we get the following results by RK45 
method: 



e 


N 


a 


\x"(0) -a\ 


\l-x'(T)\ 


x(T) 


1.0e-08 


189 


0.37939^ 


5164451122 


2.47e-08 


3.5e-08 


2.673055448977 


1.0e-09 


316 


0.37939^ 


5187063634 


2.04e-09 


2.9e-09 


2.673055570853 


1.0e-10 


549 


0.37939^ 


5189005442 


1.03e-10 


1.4e-10 


2.673055581319 


1.0e-ll 


961 


0.37939^ 


5189086642 


2.20e-ll 


3.1e-ll 


2.673055581757 


1.0e-12 


1688 


0.37939^ 


5189106905 


1.69e-12 


2.3e-12 


2.673055581866 


1.0e-13 


2827 


0.37939^ 


5189108438 


1.62e-13 


1.9e-13 


2.673055581874 


1.0e-14 


4634 


0.37939^ 


5189108571 


2.91e-14 


0.0e-14 


2.673055581875 



Remarks - as above. Here x"{T) = 9.03e - 18. The value of a = 0.3793981891086 - 
here, all digits are true. 

Numerical experiment for the equation with p = 0.1 c = 1/2, = 1; taking 
e = 10~ 14 we get T = 50. The proof of the existence and uniqueness result for p < 1 fails, 
since we cannot claim that the initial value problem (11.3D has a unique solution and that 
it depends continuously on a. Hence, function h can be multivalued. If one will prove the 
uniqueness, then, due to [7J p. 172, h will be continuous and all results of this paper will 
be true also for p < 1. The stability of numerical experiments cited below suggests it is 
the fact. 
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e 


N 


a 


\x"(0) - a\ 


\l-x'(T)\ 


x(T) 


1.0e-08 


142 


0.443643205985427 


2.16e-07 


4.5e-07 


48.05426086324 


1.0e-09 


256 


0.443643403844908 


1.78e-08 


3.7e-08 


48.05428058120 


l.Oe-10 


466 


0.443643420192529 


1.49e-09 


3.1e-09 


48.05428221034 


l.Oe-11 


839 


0.443643421402885 


2.80e-10 


5.8e-10 


48.05428233096 


1.0e-12 


1505 


0.443643421660499 


2.25e-ll 


4.7e-ll 


48.05428235664 


1.0e-13 


2669 


0.443643421681506 


1.49e-12 


3.6e-12 


48.05428235873 


1.0e-14 


4922 


0.443643421683245 


2.45e-13 


2.2e-16 


48.05428235890 



Remarks as above. Here x"(T) = 1.02e - 15. The value of a = 0.443643421683 - 
here, all digits are true. 



5 Conclusions 

The authors know that our computation of the value of the second derivative of the 
solution are not more exact than others. However, the proposed method gives a possibility 
of controlling errors and it is very simple. We hope a similar approach can be applied for 
more general equations as x'" + f(x) • g{x") = with qualitative assumptions on functions 
/ and g. 
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